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Abstract 

High momentum transfer electrodisintegration of polarized and unpolarized 
deuterium targets, d(e, e'p)n is studied. We show that the importance of final 
state interactions FSI, occuring when a knocked out nucleon interacts with the 
other nucleon, depends strongly on the momentum p n of the spectator nucleon. In 
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particular, these FSI occur when the essential contributions to the scattering am- 
plitude arise from internucleon distances ~ 1.5 fm. But the absorption of the high 
momentum 7* may produce a point like configuration, which evolves with time. In 
this case, the final state interactions probe the point like configuration at the early 
stage of its evolution. If the point like configuration is still small after propagating 
about 1.5 fm, the FSI are suppressed. The result is that significant color trans- 
parency effects, which can either enhance or suppress computed cross sections, are 
predicted to occur for ~ AGeV 2 > Q 2 < 10 (GeV/c) 2 . We suggest searching for 
color transparency phenomenon by examining ratios of experimentally measured 
quantities. Possible theoretical uncertainties of the calculations, including those due 
to the deuteron wave function and relativistic effects, are found to be small. 

PACS: 25.30.-c; 25.30.Fj; 25.45.-z 

1 Introduction 

Color transparency (CT) and color coherent effects have been recently under intense 
experimental and theoretical investigation. The (p,pp) experiment of Carroll et al.[|l| 
found indications of color transparency while the NE18 (e,e'p) experiment ]2[ found 
no such indications. The appearance of color transparency depends on formation of a 
point-like configuration (PLC) by hard scattering. The Q 2 of the NE18 experiment (1 
< Q 2 < 7GeV 2 ) seem to be large enough to form a small color singlet object. The reason 
why significant color transparency is not observed is the rapid expansion of PLC to nearly 
normal size (a nearly normal absorption) at the relatively low momenta of the ejected 
protons|| ||. Thus models of color transparency which reproduce the (p,2p) data and 
include expansion effects predicted small CT effects for the NE-18 kinematics, consistent 
with their findings, see the discussion in Ref. 0. 

Therefore a new strategy of using light nuclei and double scattering as a means 
to investigate whether small objects are produced at intermediate Q 2 was suggested 
in Ref. |]]. The idea was that using the lightest nuclei would allow the suppression of 
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expansion effects, while employing kinematical cuts to suppress the plane wave Born 
contributions and to enhance the contribution of double scattering interactions would 
increase the deviations from the Glauber model. That work concerned He nuclei target, 
in the present we use an improved formalism to investigate the use of a deuteron target. 

The deuteron is the best understood nuclear system, with a wave function determined 
experimentally in a wide momentum range @. Therefore finding a method to use the 
deuteron to investigate color coherent -color transparency effects is both challenging and 
potentially rewarding. The average separation between the nucleons in the deuteron is 
large, so the deviations from the plane wave impulse approximation (PWIA) computed 
within the ordinary Glauber approximation are usually very small, much smaller than 
for a He target. 

However a recent analysis of hard quasielastic d(p, 2p) reactions within the framework 
of conventional Glauber approximation || has found substantial Glauber screening effects 
for transverse spectator momenta less than 300 MeV/ c and small longitudinal spectator 
momentum. In particular, the interference between the Born term and rescattering 
amplitude in which the knocked out nucleon interacts with the spectator nucleon, causes a 
reduction in the cross section. This is called a screening effect. Including the square of the 
rescattering amplitude, which we call the double scattering term, increases the computed 
cross section. Color transparency leads to a suppression of the rescattering amplitude, 
predicting specific upward and downward changes in the cross section, depending on the 
kinematics which control the relative sizes of the screening and double scattering terms. 
Thus observing such effects which occur for different values of the spectator momenta 
seems to be an effective method to search for color coherent effects using the deuteron. 

Substantial screening has been observed also in the recent calculation || of the pro- 
cesses d(e, e'p) within the nonrelativistic Glauber approximation. The present paper uses 
formulae of the Glauber approximation deduced from the Feynman diagrams. As a re- 
sult the predicted angular dependence of the cross section at higher spectator momenta 
is noticeably different from that in Ref. 0. In particular, the position in the spectator 
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angular distribution, corresponding to the maximal final state interaction depends on 
the momentum of the spectator nucleon. At Q 2 ~ l(GeV/c) 2 our calculations agree 



well with calculations of Ref. [10], |TTJ in the kinematics where the double scattering term 



dominates. According to Ref. [10] in this kinematics the contribution of meson exchange 
currents and isobars in the intermediate is small. This is another example that the 
Glauber approximation is a good starting point in the hunt for the CT phenomenon. 

We stress that exclusive processes of high momentum transfer electrodisintegration 
of the polarized and unpolarized deuteron - d(e, e'p)n have an important advantage. 
One may choose special kinematical conditions for the spectator nucleon to control the 
essential internucleon distances (~ 1.5 fm) for final state interactions (FSI). Therefore it 
is possible to probe the PLC at the early stage of its evolution, where expansion effects 
are minimal. The chosen kinematics also allow the separate investigation of the different 
implications of CT for the screening and double scattering terms. Another limitation 
on the kinematics is the requirement of reliability of the nonrelativistic description of 
deuteron and smallness of competing nuclear effects which might diminish CT effects. 
The major price for fixing kinematics is the drop in the cross section, which should not 
cause a problem for high intensity accelerators. 

The paper is organized as follows. In section 2 we calculate the cross section of the 
d(e, e'p)n reaction within the framework of the conventional Glauber approximation [|T^ . 
Polarized and unpolarized targets are studied. 

In section 3 color transparency effects in the d(e, e'p)n reaction are calculated within 



the framework of quantum diffusion model|3|] and three resonance model[13|. The three 
resonance model is solved in two ways: the Green function method of Ref. || and the 
modified Glauber approximation method, used previously to calculate nuclear vector 



meson production (see Refs. fl~4], |I5| , [16| , p"?!). 

Section 4 contains the numerical results. We first use the Glauber approximation to 
investigate the kinematic requirements necessary to ensure the most favorable conditions 
for studying color coherent effects. This analysis allows us to introduce new experimental 
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observables which are free from theoretical normalizations and are more sensitive to the 
effects of final state interactions. Using the models, discussed in section 3, we calculate 
color transparency effects for the kinematic range accessible to the present generation of 
electron accelerators. 

Section 5 considers theoretical uncertainties and the reliability of obtained results. 

The main results of present paper are summarized in Section 6. 



We consider the d(e, e'p)n reaction for kinematics in which the proton carries almost all 
of the momentum of of the virtual photon and the neutron momentum p n is in the range 
200-400 MeV I c where deuteron wave functions have been determined experimentally. [] 
The amplitude for d(e, e'p)n processes can be expressed as the matrix element of the 
scattering operator - T between the deuteron wave function and the wave function of the 
final state consisting of two outgoing nucleons: 



Here | d, s) is the position space deuteron wave function, p p , f p , s p and p n , r n , s n are the 
momenta, coordinate and spin of proton and neutron in the final state. 

Within the conventional assumptions on factorization of hard - electromagnetic and 
soft - FSI amplitudes i.e. within the Glauber approximation the nuclear scattering op- 
erator has the coordinate space form: 



where Tf[ x {Q 2 ) - is the one-body electromagnetic current operator, q is the space com- 
ponent of the photon four-momentum which we choose in the z-direction. We note that 
1 Obviously, kinematics in which the neutron is produced forward and proton sidewise is as good for 
our purposes. In fact, such a setup may have certain experimental advantages since in this case one 
detects neutron with the same momentum for different Q 2 . We thank B.Mecking for emphasizing this 
point. 



2 Glauber Approximation (GA) 





T(q,r p ,r n ) = T s (p p ,p n ,r p ,r n ) -T^iQ 2 ) ■ e 



(2) 
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factorization is violated when off-shell effects are taken into account. So we will focus on 
CT effects at small spectator nucleon momenta or in the rescattering region where only 
small momenta in the loop integral over internal nucleon momenta are important. We 
neglect the contributions of meson currents like that at Fig.fc since such contributions 



rapidly decrease with Q 2 for x ~ 1 (cf. discussion in |18| and sec. 5.5 bellow). T s (p p ,p n ) 
describes soft FSI of knocked-out proton with spectator neutron. We are studying situ- 
ations in which the momentum of knocked - out proton is larger than 1 GeV/ c. So Tg 



can be calculated within the Glauber approximation WA IK 



T s (p P ,Pm r p , r n ) = 1 - T N (b p - b n ) -<d(z n - z p ) ■ e -*A°fc>-*0, (3) 

where r = f(z, b). Here the z direction is defined by the direction of the struck proton 
momentum p p . However within the considered kinematics, where proton carries almost 
all momentum of virtual photon, we choose for z the direction of q. The profile function 
T N (b) is expressed via the nucleon -nucleon scattering amplitude (/ ) as follows: 

r N (b) = I /exp (ik t ■ b) ■ f NN {kt)^ (4) 

where the NN scattering amplitude normalized as Imf NN (k t = 0) = Otot- The addi- 
tional factor e -* A °( z P-zn) j n eq.fl^), where 

A = (E n -m) R — , (5) 

is a kinematical factor. In eq.fl^) E n - is the spectator energy, m is the nucleon mass, 
ma - is the mass of deuteron and go is the transfered energy. The factor A is obtained 
by evaluating Feynman diagrams (cf. section 5.4, eq.fl5"8D) and accounts for the nucleon 
recoil usually neglected within the conventional non-relativistic Glauber approximation. 
This factor accounts for the fact that the variable fc °~ fcz but not k z (fig. lb) is conserved 
in two-body high energy collisions (cf. f2T||). 



2 To simplify the discussion we neglect small corrections due to the charge exchange process np —>■ pn. 
Since the charge exchange amplitude is predominantly real and of spin flip character (which itself is 
small), it practically does not interfere with the main amplitude. 



Inserting eqs.(|3]) and ([|) into eq.(0), we obtain the amplitude for the d(e, e'pjn reaction 

as: 

M = T e H m (Q) f fd 3 r p d 3 r n e^(r p ,s p ,r n ,s n \d,s) 



1 - T(b p - b n )e- lA °^-^&(z n - z, 



TfTiQ) J J d\d 3 r n e- i ^e- i ^(f P7 s p ;r ni s n \ d, s) 



x hl e *' ( ^ n) ■ /^A)^ 6 "^ ^ 6 ^ - *)> ( 6 ) 

where we define the momentum of initial proton in the deuteron state as Pi = p P — <?• 

To separate center-of-mass motion we introduce the center-of-mass (R cm ) and relative 
distance (r) coordinates: 

T*p Rem ~)~ 

^ = ^--^ ( 7 ) 

Using Eq.(^) allows us to simplify Eq.@ as: 

M = Tp m (Q) J d 3 R cm e lli --^-^-^ || t Pr(? p ,s p -,7 f n ,s n \d,3)e-fr*efc* 



d 3 r&(-z)(f p , s p , r n , s„ | d, S) • e-^-e^I / e *-S-«A°« . fN^fh 



2iJ v y (2 



7T 
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- ^ // ^e(-z)(r P) r n; a „ I d, s) ■ e-m.-m-iiP^Aoy^NNfy^^ (g) 

The integration over i? cm ensures the momentum conservation: p n = —pi = q ' — p p . In 
the last part of eq.(§) we omit the factor (27r) 3 5(pd — p p — p n ), which will be included 
in the definition of cross section. (To simplify formulae we use the Laboratory frame, in 
which the momentum of of the deuteron, is zero.) 

The integration over the relative coordinate r leads to the result: 

M = (2rr)i • T; m (Q) U Pi , s p , s n \ drf-^-J ( Pi , s p , s n \ d, ^f NN {k t )^^ , (9) 



where p\ = (pi Z + A°,p it — k t ) and (p[, s p , s n \ d,s) - is the deuteron wave function 
in momentum space representation. The additional factor | in the second term is a 
consequence of Q(—z) functions in eq.(§). (The accuracy of last replacement is discussed 
in sec. 5. 4). 

Taking the square of the modulus of the amplitude and summing over the proton and 
neutron polarizations we obtain: 

|| M || 2 = MM] 



(2vr) 3 • |F— (g 2 )| 2 ■ I pjfoft) - Tie- J ^pM, Pi )f Nm (k) 

Hi,Pi2)f NN (hi)f NN *(h 2 )\, (10) 



16 7 J {2tx) a 

where p~l x = (p i2 + A°,p it - k tl ), p i2 = (p iz + A ,^ - fc t2 )- 
In eq.fllOD we introduced the deuteron density function: 

p s (k ll k 2 )= J2 ( k 2,s P ,Sn \d,s){d,s\ ki,s' p ,s' n ). (11) 

Sp,s' p ,Sn,s' n 

In this Glauber or distorted wave impulse approximation (DWIA) the d(e, e'p)n cross 
section can be expressed as follows: 

dE e Xd* Pp = ' D *&P»P«) ■ S ^ ~ Md - E *- E ^ ^ 

Here a en is the cross section of the electron scattering off a bound proton (up to the 
flux and proton recoil factor) f\ and the decay function Dd(q,p p ,p n ) represents the joint 
probability for the initial proton in the deuteron having Fermi momentum - and for 
the final state having a proton and neutron with a momentum p p and p n (for more details 

see 



22 



3 In principle, one should use the light-cone quantum mechanics of the deuteron plj | to calculate the 
cross section. However for the effects dominated by the contribution of small nucleon momenta in the 
deuteron of interest here, the difference between predictions of light-cone and nonrclativistic formalisms 
is small (see Sec. 5.1). 
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If the deuteron is polarized, the cross section depends on D d (q,p p ,p n ) with 
D° d {q, p p , p n ) = P g d {pu Pi) ~ Tle^. J p g d { Pi , p[) ■ f NN {k 



d 2 h 



(2tt) 2 

/ ^(Ri.Po) " f NN (ka) ■ f NN *(k t2 ) ■ i k ti k l . (13) 



16r BVU,rM/ J v "" iy J V " M/ (27r) 2 (27r) 2 ' 
The polarized density matrices of deuteron can be expressed via the s and d - wave 



components introducing the polarization vector a according to Ref . P5 

u(ki)w(k 2 



Pd(ki,k 2 ) = u(k 1 )u(k 2 ) + 



3\ko ■ a\ 2 



A 9 



>/2 



+ 



3\ki ■ a\ 



u(k 2 )w(ki 



+ 



'9 (A;i • a)(k 2 ■ a)*(ki ■ k 2 ) 3\ki-a\ 



k 2 k 2 



2 fc? 



(14) 



where the components of a are defined through the deuteron spin wave function as: 



^ 10 = z • a 2 , ip 11 = i=(a x + ia y ), ip L 1 = —j={a x — ia 



.ii 



y)i 



(15) 



where tp 1 ^ is the projection of the deuteron's spin on the the \x direction. The unpolarized 
deuteron density matrix follows from Eq.flTj]) as: 

p d (h, fe) = | £ P8(*i, h) = u{k x )u{k 2 ) + w{k x )w{k 2 ) ■ {^^f- ~ ~j ■ (16) 

3 Coherent Effects in <i(e, e'p)n Scattering 

Theoretical analysis shows that in realistic models the absorption of a hard photon leads 
to the formation of a point like configuration, which undergoes a reduced interaction 
with other hadrons, because of its small size and its color neutrality [p^j . To estimate 
the expected effects of color coherence we consider two different models which account 
for the formation of the PLC and their evolution to the normal hadronic state: quantum 



diffusion model of Ref. and the three state model of Ref. |13 
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3.1 Quantum Diffusion model (QDM) 

The reduced interaction between the PLC and the spectator neutron can be accounted 
for by introducing the dependence of the scattering amplitude on the transverse size of 
PLC. However we consider energies that are far from asymptotic, so the expansion of PLC 
should be important. This feature is included by allowing the rescattering amplitude to 
depend on the distance from the photon absorption point. Including these effects leads 
to the modified deuteron decay function D d (q,p p ,p n ): 

Di(q } p P ,Pn) = p s d {u(p i ),w{p i ),u{p i ),w(p i )) 

-Tie- [ pi{u{p t )MPi)Mp'yMp\Y)f NN {kt ,H " 



I J (27T, 



2 



+7 1 ^(^(^)^fe2),u(K)>(K))|%r 7v (fc tl )r Ar *(^)^ // 



(2tt) 2J v lLJJ v " y (2vr) 2 ' 

(17) 



with the same form of p s d as in Eq.flHD but with modified radial wave functions: 

1 r fPLC,N( u f)2\ 

m = / <r)Y (r) • ; f k ' Q } e{-z)e-«*<Pr 

(2ll)2 J J [Kt) 

1 r fPLC,N( u f)2\ 

m = -At / w(r)Y 20 (r) ■ 1 J™ V > Q(-z)e**<Pr. 



(27T)* J " ' " ' f NN (kt) 

Eq.dTTD is derived in the same way as the Glauber formulae in the previous section. 
The only difference is to include the coordinate dependence of the rescattering amplitude. 
To calculate the deuteron decay function in the quantum diffusion model we model the 
amplitude of the PLC — N scattering in a form|| consistent with the optical theorem: 

f PLC > N (z, k t , Q 2 ) = ia tot (z, Q 2 ) ■ el* • ^t^lpI^A , (19) 

where b/2 is the slope of elastic iVTV amplitude, Gjv(t) (~ (1 — t/0.71) 2 ) is the Sachs 
form factor and t = —k 2 . The last factor in eq.(|l9D accounts for the difference between 
elastic scattering of PLC and average configurations, using the observation that the t 
dependence of da h+N ~' h+N / dt is roughly that of ~ G 2 (t) ■ G%(t) for not very large values 
of t and that G 2 h (t) w exp(R 2 h t/3). 
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In Eq. (PUD <Jtot{li Q 2 ) is the effective total cross section of the interaction of the PLC 
at the distance I from the interaction point. The quantum diffusion model |§ corresponds 
to: 

v t ot(i,Q 2 ) = cjtot | (j- h + (r * |^ )2) (1 - jjj 0(l h -l) + 0(1 - / h )J , (20) 

where l h = 2p f /A M 2 , with A M 2 = 0.7 - 1.1 GeV 2 . Here (r t {Q 2 ) 2 ) is the average trans- 
verse size squared of the configuration produced at the interaction point. In several 



realistic models considered in Ref. [24| it can be approximated as ^ ~ 1 G ^ 2 for 



Q 2 > 1.5 GeV 2 . Note that due to effects of expansion the results of calculations are 
rather insensitive to the value of this ratio whenever it is much less than unity. 

3.2 Three-State Model 
3.2.1 Summary of the model 

To evaluate matrix element of operators Tg and T# in eq.(^) we use the three-state model 



of Ref. |13| . The basic assumption of this model is that the hard scattering operator T# 
acts on a nucleon to produce a non-interacting point-like configuration \PLC) which is 
a superposition of three baryonic states: 

T H \N) = \PLC)= F m , N (Q 2 )\m), (21) 

m=N,N* ,N** 

where F m ^(Q 2 ) are elastic (m = N) and inelastic transition form factors in this model. 
It is important to notice that the transition form factors used here cannot be taken 
directly from data. This is because states iV* and iV** are effective states which, in fact, 
represent a number of actual physical states. In the following analysis we assume for 
certainty that all form factors have the same Q 2 — dependence. We neglect also possible 
spin effects in the form factors and in the operator X#. 

Color transparency is introduced in this model as the condition of lack of FSI at the 
point where PLC is produced: 

T S \PLC) = 0. (22) 
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It is this condition which distinguishes models with color transparency from models which 



include production of resonances and/or continuum in intermediate states ||25||. Within 
the approximations discussed above T5 is equal to a matrix U which is the most general 
3x3 (Hermitian) matrix that annihilates the PLC: 



U 



a 



tot 



( 



\ 



Fn*,n 



\F N ,N\ 2 -e*F N , N F^^ N -tM\F N ^ N \ 2 

F N* tN F N**,N 
y\F N *, N \ 2 -\F NtN \ 2 +2Re(e*F N , N F^ t N ) 

\Fn**,n\ 2 '~ 



(23) 



Taking account of the data on pp and pd diffractive scattering helps to restrict further 
the parameters of this matrix|T3|. We use the matrix U from that work. 



3.2.2 Green Function Method 

In this section we use the Green function method developed in Ref.0 to calculate CT 
effects in the double scattering reaction off 3 He. Within this method, the matrix element 
for the process d(e, e'p)n which includes final-state interactions, is equal to 



M s {p p , s p ;p n , s n ) = (p p , s p ;p n , s n \T H + T s G T H \d, s), 



(24) 



where p p ( n ) is the momentum of the recoil proton (neutron), s p ( n ) is the projection of 
the recoil proton (neutron) spin and s is the projection of the deuteron spin. Th is the 
hard scattering operator - the electromagnetic current. Go is the free propagator, and Tg 
accounts for the soft final-state interaction. The matrix element M. s is also a function of 
the four-momentum transfer Q 2 and the Bjorken variable x (or equivalently the virtual 
photon three-momentum q and its energy q°). However, this dependence will be omitted 
in the formulae to simplify notations. It is convenient to define 



■M-oiPpi S pj Vni s n) — (Pp: s p] Pm s n\Tf[\d, S) , 

M{(p p , s p ] p n , s n ) = (pp,Sp]p n ,s n \T s GoT H \d,s). 
The matrix element for the hard scattering operator is given by (cf eq.(pl~P) 

(m,f p ,s p ;n,f n ,s n \T H \d, s) = F mtN (Q 2 )e ig ' Tp (r p , s p ;r n , s n \ d,s). 



(25) 
(26) 

(27) 
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Here m labels components in the PLC (N, N* , N**). The Green's function operator is 
diagonal in the hadronic mass eigenstate basis : 

e jpm|r>-r^| 

(m,r p ,s p ;n,r n ,s n \G\m,r p ',s p ;n,f n ',s n ) = -— — —r5 {3 >(f n - r n '). (28) 

4ir\r p — r p '\ 

Here p m is the momentum of the m'th component of the proton wavepacket. This is 
given by: 

Pmz Q Vnz 

Q\l-\) + Ml-Ml 

Vnz ~ 2g ■ W 

Here p' nz is the momentum of the neutron spectator in the intermediate state (cf. discus- 
sion in sec. 5.4). The above expression is obtained using energy-momentum conservation 
in intermediate states, valid in the semiclassical approximation used here. The eq.([48|) 
shows that when \PLC) state is modeled as a superposition with different masses, the 
longitudinal momentum of nucleons in the intermediate state differs from the longitudi- 
nal momentum of registered nucleons. This is another important difference of our three 
resonance model [DJ from the model discussed in Ref.|25[ . 



Here z is the direction of the 7* momentum, x is the Bjorken variable, and we take 
for simplicity M p = M n = and neglect the binding energy of the deuteron. The 
quantity is the squared mass of the m'th component of the wavepacket state. 

The matrix element of the final-state interaction operator T5 is parameterized ac- 
cording to eq. (p3|) with 



(rp, Sp, r n , Sj-jT^glTTi, v py Sp, m , r n , s^) U m ^ p 5^ (Tp fn)- (^0) 

For the simplicity we neglect here the dependence of scattering amplitudes on momentum 
transfer (See however discussion in sec. 3.2.3). 

We put all of the pieces together, insert complete sets of states, perform the integrals 
over the delta functions, and arrive at the following expression for the matrix elements: 

M s Q (p p ,s p ;p n ,s n ) = F PiP (Q 2 ) Jd%d 3 r n e- l ^e- l ^"e^(r' p ,s 
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Ml{p p ,s p ;p n ,s n ) = i\p p \J2U P>m F m , p (Q 2 ) I '^re-'W^ 

m 

xe p - ^rjyvi s v'i r nS n d,s), (32) 

47r|r p -r n | 

We introduce the relative and center-of-mass coordinates according to eq.(^), with this 
transformation, d 3 r p d 3 r n = d 3 R cm d 3 r. The integral over d 3 R cm leads to momentum 
conservation: q = p p + p n in the Lab. frame of the deuteron. The remaining matrix 
elements are expressed purely in terms of the relative coordinate, r. 

The non-relativistic deuteron wavefunction |7J can be written in terms of the s- and 
d-state functions u and w as, 



(r | d,s) 

where the tensor operator 



u(r) w(r) 

7^012 (r) 

r y/8r 



\s), (33) 



c ,^ 3a p -ra n -r _ _ 

S12 (r) = -3 <V<Xn- (34) 

As a result of the angular integration in d 3 r (in particular, the integration over az- 
imuthal angle ) s p + s n = s. It is, thus, convenient to define the quantities 

M s m (p p ,p n ) = M s 0{1) (p p ,s p ;p n ,s n ). (35) 

When s = ±1 then s p = s n = ±|. When s = 0, then .Mom is symmetric combination 
of s p = — s n = I and s p = —s n — —\- The amplitudes for the production of final-state 



proton are 



M s =±l (p p ,p n ) = 2irF pp (Q 2 ) [Uo + Wo] , (36) 

M s =°(p P ,Pn) = 2-kF pp (Q 2 )[U -2W ], (37) 

M{ =± \p p J n ) = 2TrJ2U p , m F m , p (Q 2 )[U m + W m }, (38) 

m 

M s r°(p P ,Pn) = 2nJ2U P , m F m , p (Q 2 )[U m -2W m ], (39) 
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where 



Wo 

Urn 



sin p n r 



2 / drruir) 
Jo p n r 



V2 / dr rw(r) 
Jo 



sm p n r cos p n r 

T O- 



, smp n r 

^3 rr* 3 

fn 



dr e ippmT u(r) mg , 
47r Jo gr 



\/2 f 00 

/ dre Wpmr w(r) 

An Jo K ' 



sm gr cos gr sm gr 



qr 



iy 2 



^3 ^ 3 



(40) 
(41) 
(42) 
(43) 



and where g = |g| and p n = \p n \. Note that Wo = for p n = 0. At very high energies, 
p Pm ?z p p is independent of m, so that U m and W m become independent of m. In this 
limit the matrix element J\A.\ is proportional to J2 m U Pim F myP = 0. This is the limit of 
full color transparency. 

The formulae of the conventional Glauber approximation can be obtained from the 
above equations by retaining only the proton component in this sum. 



3.2.3 Modified Glauber Approximation Approach 

In this subsection we solve the three-state model using the method developed by Yennie 



14}| , Kolbig and Margolis |LJ| and Bauer [1J|. The major difference from the standard 
Glauber approximation is the modification of the profile function by substituting /nn—* 
f nm , where f nm is the m-particle production amplitude and the phase factor due to the 
difference of the masses of particles n and m|l6|]. This differs from the previous subsection 
and provides an opportunity to include a non-zero range of the interaction (the slope b m 
of the amplitude f nm ). 

Thus we can use equations of section 2, and modify the profile function in eq.@ to 
calculate the d(e, e'p)n processes within the Three State Model: 



T N (b p -b n 



E 

m=N,N*N* 



j r Fm,p %M ikf(r v -r n ) -iA m , v (z v ~z n ) ffij 

P ' m F n J (2tt) 2: 



where 



A 



m,p 



. ]\,j 2 — M 2 

\p m \ - \p p \ = A -^M 2 N +p 2 p -M 2 m -\p p \ = A - " — ■ 



(44) 



(45) 
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In eq.(fH|) b m is the slope factor of the transition amplitude. For the present analysis 
we will assume that slopes of all transition amplitudes are the same as for NN scattering. 
The factor e lA m, P (z P -z n ) acc0 unts for the phase shift due to the different masses of PLC 
components. The additional factor accounts for the different coefficients for \PLC > 
wave components. The main difference of eq. fl44j) from the formulae of the Generalized 
Glauber approximation |L4], [15| , |16[ |P7|| is that eq.( [i4l) should satisfy the sum rules eq.(^) 
in the point where PLC is produced (z p — z n = 0) and when t — * 0. Besides the effects 
of the m-dependence of t m i n are effectively included in the sum rule eq.(|2"2|). 

Using modified profile functions in the formulae for the transition operator in eq.(fj) 
we obtain: 

M = F P:P (Q 2 ) J d 3 r(f, s p , s n | d, s) {e-** 

-VXI Ph U P^e^Q(-z) ■ e -^ e -^^A . (46) 
^ m J y^i ^P,P ) 

The ^-component for the momentum of knocked out proton in intermediate state is : 

V Z = p iz + A m>p , (47) 
Finally we obtain the following expressions for the deuteron decay function: 

D d (q,p P ,p n ) = ffoiM) ~\ E / Afatf) " U p , m ^ -^-^k 

Z m=P,N* ,N** J P P>P V Z7F / 

^h 1 ^ p ' m ' fl {2,f{2,r (48) 

where pf' = PT' (p™,Pu ~ h), Pa = Pa (p™, Pit ~ hi), Pa = PS(p™,P;t ~ hi), and 
density matrices p s d is same as defined in eq.([TJ|). 



4 Numerical Results 



To estimate effects of final state interactions in the d(e, e'p)n reaction we compare the 
predictions of different models (considered in sections 2 and 3) for the differential cross 
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section and the tensor polarization asymmetry with predictions of the plane wave Born 
impulse approximation. 

The manifestations of color transparency depend strongly on the relative importance 
of the screening (interference term between Born (Fig. la) and FSI amplitude (Fig. lb)) 
and double scattering terms (square of FSI amplitude). To separate different implications 
of CT it is necessary to distinguish kinematical conditions where one of these terms is 
dominant. Thus we analyze the deuteron decay function in eq.([H|) defined in Glauber 
approximation, for the separate cases of unpolarized and polarized deuteron targets. 



4.1 (e,e'p) scattering on unpolarized deuteron target 

Consider first the unpolarized deuteron case. We define the transparency T as the ratio 
of the measured cross section (or calculated cross section with FSI) to the one calculated 
in the PWIA |], |, |, |, |§: 



T(0 2 v v ) = ad(e ' £/p)ra(g ,Pp,Pn) (49) 

a d(e.e'pn)W iPpiPn) 

Within the Glauber approximation, the cross sections are controlled by the decay function 
in eq.(Ji"3|) with the unpolarized density matrix of eq.(|T^). In Fig. 2 the dependence of the 
transparency - T(Q 2 ,p p ,p n ) on the spectator neutron angle 6 n is presented for different 
values of the spectator momenta - p n at fixed Q 2 = 6 (GeV/c) 2 . The figure demonstrates 
that the final state interaction is maximal for 9 n « 90° at smaller spectator momenta 
(< 200 MeV/c) and the position of maximal FSI shifts to slightly lower spectator angles 
with increase of p n . The shift is the consequence of kinematic factor A of eq.((5|) which 
accounts for the fact that for rescattering of energetic particles the E n —p z n component is 



conserved rather than p z n (see eq.(^) in section 5.4 below). As follows from eq. (|58|) the 



maximal FSI effect occurs at x = = 1, which provides the zero value of longitudinal 
momenta for intermediate spectator. Below we will refer to the exclusive kinematics 
corresponding to x = 1 as the perpendicular kinematics. 

The fig. 2 shows also that the final state interaction contributes differently to the re- 
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suiting cross section at different spectator momenta. At p n < 200 MeV/c (Fig2.a,b) FSI 
is dominated by the screening effect (the second term in eq.fljjD). At p n > 300 MeV/c 
(Fig.2.c.d) the double scattering term (the third term in eq. ([T3|) ) plays the dominant role 
at the x ~ 1 and it tends to diminish the screening effects (dash-dotted curves in Fig. 2). 
The sharp increase of the FSI with the increase of spectator transverse momenta is the 
consequence of the decrease of the values of actual impact parameters contributing to 
NN rescattering. Note that in calculations within the Glauber approximation the pn 
scattering amplitude have been estimated using the relation: 

f pn = • (i + «n)e 6 " /2 * (50) 
where ot n = Ref /Imf and parameters used are those of Refs.|26], |27|, ^8], ^9 . 



The above analysis suggests that color transparency effects, which reduce the FSI 
(rescattering amplitude), will have different consequences for different intervals of spec- 
tator angles and momenta. 

Results presented in Fig. 3a indicate that transparency for the scattering off the 
deuteron in perpendicular kinematics (x = 1) should change quite dramatically with 
p n from T{p n = 0) ~ 0.97 down to T(p n = 0.2 GeV/c) ~ 0.5. The decrease of T shows 
that FSI, in this kinematical range, is mostly a screening effect (see Fig. 2). Note that 
current NE-18 results@ for T d were obtained under assumption that T does not depend 
on p t in the studied momentum range p t < 300 MeV/ c. Therefore, to compare our results 
with these data it is necessary to perform a new analysis of the data with a more realistic 
nuclear model where T(p t ) decreases with p t . 

However at p n > 300 MeV/c, x — 1 the transparency T starts to increase with 
spectator momentum (Fig.3b), since in this kinematics (see Fig.2) the double scattering 
term dominates. 

As a consequence of different role of FSI in the considered two kinematics, CT effects 
will have nontrivial manifestations. CT will increase T as compared to GA predictions 
at p n < 200 MeV/c, x — 1 (see Fig. 3a- curves "I, II, III") and will decrease T as compared 
to GA predictions at p n > 300 MeV/c, x = 1 (see Fig.3b- curves "IV, V"). 
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The different pattern of CT effects for these two regions shows that the ratio of the 
cross sections: 



, a(p nl « 300, 400Me1//c) 



R{Q ^ Pn2) = a( Pn2 ~ 200MeV/c) (51) 
should be rather more sensitive to CT phenomenon. This quantity is more convenient to 
search for CT since it represents the ratio of directly measured experimental quantities, 
and does not require additional normalization to the corresponding PWIA calculation 
(as in eq (|49"1) ) . In Fig.4 we present curves for the same kinematics as in Fig.3 where 
calculations accounting for CT effects are normalized to the corresponding Glauber ap- 
proximation calculations. We calculate also the Q 2 dependence of the ratio- R (curves 
labelled by the boxes). One can see that CT can modify R by as much as 30% for Q 2 as 
low as 6 — 10 (GeV/c) 2 . Note that this occurs in the kinematical region where competing 
nuclear effects in the deuteron are small and are under good theoretical control (see sec. 
5). 

In Fig. 5 we compare the predictions of different CT approximations for transparency 



defined as in eq.(p9|). Comparison of calculations within the three resonance model 
shows that the Green function formalism (see sec. 3. 2. 2) predicts smaller DWIA cross 
section(Fig.5 - solid lines labeled by boxes) than the modified GA approach (see sec. 3. 2. 3). 
This difference is caused by the slope of the NN scattering amplitude in the modified 
GA approach. Including this effect allows the transverse internal nucleon momenta in 
the deuteron to be smaller in the modified GA approach. This is a qualitatively new 
feature. However, the effects of the finite range of interaction do not exceed 10% in 
the discussed kinematical range p n < 300MeV/c. In the case of CT, for spectator 
momenta (p n < 200MeV/c), neglecting the slope of the soft rescattering amplitude 
predicts approximately the same effects as modified GA (Fig. 5a - dashed curves labeled by 
boxes). For larger spectator momenta (p n > 300MeV/c) neglect of slope for rescattering 
amplitude enhances CT effect (Fig. 5b - dashed curves labeled by boxes). This is because 
for higher spectator momenta rescattering occurs effectively at smaller distances than 
in the case of the modified GA. These considerations demonstrate that up to deuteron 
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internal nucleon momenta ~ 200MeV/ c one can neglect the finite sizes of the rescattered 
nucleons, while at higher momenta when internuclear distances become comparable with 
nucleon size one needs to account for the form factors of the interacting nucleons. 

The energy dependence of final state interactions of the d(e, e'np) reaction involving 
a spinless deuteron target has been recently discussed in Ref. f32|. That paper studies 
ratios of cross section integrated over for various values of Bjorken x and uses the triple 
pomeron mechanism of diffraction, which is not applicable for the energies considered in 
this paper. 



4.2 (e,e'p) scattering off the polarized deuteron 



The possibility to employ a polarized deuteron target to investigate color coherent effects 
is quite tempting. Use of different polarization states enhances the role of the ci-wave 
component of the deuteron wave function and so that smaller space-time intervals are 
probed. This results in tagging of the PLC at the early stage of its evolution to a normal 
hadron. 

For numerical estimates we consider the tensor polarization T 2 q measurable in elec- 
trodisintegration of the polarized deuteron. We define T 20 as: 

T 20 = ^ 0(1, 1) + (7(1, -1) - 2 • <t(1, 0)) = a ep -D^(q,p p , Pn )-S{q -M d -E p -E n ), (52) 

where a(s, s z ) = dE ^"* d 3 p > s anc ^ Sz are ^ ne s P^ n an< ^ z com P onen t of the deuteron. 
The decay function is defined according to Eq.(|i"3"D with the tensor polarization density 



function as: 

pto(h,k 2 ) 



3 \P d 



7 ' V \k 1 ,k 2 )+ P p {a ^\k 1 ,k 2 )-2.p^(k 1 ,k 2 ) 



1 / -^=(a x +ia y ) 
P, 

kn. .1 u(ki)w(k 2 ) 



3 ■ -2z 

k 2 



V2 



+ 



k( 



u(k 2 )w(ki) 
V2 



+ 1 



(h ■ k 2 ) ■ ((An ■ k 2 ) - 3k lz k 2z ) *4 

ft l rv 2 



k 2 k 2 



lj w(kx)w(k 2 ). 

(53) 
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Eqs.(|T4]) and ( pT5|) are used in obtaining eq. 

For further calculations we define an asymmetry as the ratio of cross sections from a 
tensor polarized and unpolarized deuteron: 

^ Pn) = Tw/ dE e ,dn e ,d*p p ~ DT p ( q ,p P ,p n y (54) 

where the unpolarized decay function - D^ np (q,p p ,p n ) is defined according to eq. fli~3]) and 



Before considering manifestations of CT for T20 and Ad it is worthwhile to exam- 
ine some properties of the tensor polarization density matrix eq. flo"3"|) . Considering the 
scattering from unpolarized deuterons (section 4.1) demonstrates that FSI dominates for 
nearly perpendicular kinematics: p s _L q. Thus in Fig. 6 we compare P2 {k,k) at k z = 
with the corresponding unpolarized deuteron density function pd{k, k) = u(k) 2 + w(k) 2 
calculated with different (Paris and Bonn) NN potentials. Fig. 6 illuminates several re- 
markable properties of the tensor polarization density functions. First, it practically 
coincides with unpolarized density function for k « 300 MeV/c, where u(k) ~ — 
This property reflects the fact that the s- partial wave in the deuteron falls with k, and 
changes sign at k « 400 MeV/ c, while <i-wave grows with k from negative minimum at 
k fa 100 Mev/c. Thus in some range of momenta, the influence of the (i-state is compa- 
rable (or larger) than that of the s-state (for a more detailed discussion see 0). Note that 
we follow here the convention of Refs. |3(], [3T| which define the deuteron wave function so 
that w(k) < at small k. The opposite convention, corresponding to w(k)>0 is adopted 
e.g. in Refs. 0, p2|). This equality of polarized and unpolarized density functions means 



that in this kinematical range the asymmetry calculated according to eq.(p4|) will be 
close to unity, for perpendicular kinematics provided there is no final state interaction 
(i.e. in the plane wave Born approximation). This prediction is practically independent 
of the NN potential used to compute the deuteron wave function (see Fig. 6). Besides, 
we will see in sec. 5 that other effects such as the EMC effect, off-shell and relativistic 
effects are small for our kinematics. So deviations of the asymmetry from unity originate 
predominantly from the influence of final state interactions. 
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The next important feature of the tensor polarization density function is that it 
decreases with decrease of internal momenta (Fig. 6), and we find that average momenta 
in the integrals for the final state interactions in the perpendicular kinematics are 100 — 
150 MeV j c. Therefore, within the framework of the conventional eikonal approximation 
large effects of FSI are expected, since rescattering occurs at rather small internucleon 
distances ~ 1.5 fm. This dominance of FSI by small internucleon distances indicates 
that the effect of PLC evolution to a normal size nucleon will be reduced. Therefore one 
should expect larger sensitivity to CT effects. 

In Fig. 7 we present a three-dimensional plot of the asymmetry Ad as a function of the 
spectator momenta and the polar angle, calculated within the PWIA (Fig. 7a) and the 
Glauber approximation (Fig. 7b). This figure clearly demonstrates the large influence of 
final state interactions within the framework of Glauber approximation for perpendicular 
kinematics with p n w 250 - 350 MeV/c and 6 n » 60 - 90°. 

In Fig. 8 the Q 2 dependence of Ad is calculated using the PWIA, Glauber approx- 
imation and CT models considered in text, for different values of spectator transverse 
momenta. The sensitivity to CT effects is noticeable. 

In Fig. 9 the Q 2 dependence for both the ratio R{Q 2 ,p n i,Pn2) defined according to 
eq . fl5l|) and asymmetry Ad defined according to eq.flSiD calculated for different values 



of AM 2 , which determines the coherence length within the QDM prediction of color 
transparency. Here we restrict Q 2 to the range available at CEBAF in the near future 
(Q 2 ~ 6 (Gev/c) 2 )^. One can see that already in the E e = 4 - 6GeV CEBAF run, it 



is possible to obtain important constraints on the parameters of CT. 

5 Theoretical uncertainties 

In this section we consider several theoretical issues which might influence the reliability 
of our interpretation of the measured cross sections. 
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5.1 Relativist ic motion of target nucleons 



The fact that high energy processes develop along the light-cone can be taken into ac- 



count within the framework of the light-cone mechanics [21], 0- The light-cone calcu- 
lation of processes involving a deuteron target is more straightforward than that using 
nonrelativistic quantum mechanics. The parameter which characterizes the importance 
of relativistic motion of nucleons in the deuteron is the difference between the deuteron 
internal momentum defined in the laboratory reference frame (p n ) and in the light-cone 



reference frame 21 



m 2 + pl t 



— m z 



(55) 



\| a(2-a) 

where a = is the light-cone fraction of deuteron momentum carried by spectator 
nucleon. In the perpendicular kinematics of present interest a ~ 1 and it follows from 
eq . (|55"D that k n ~ p n t- Therefore the expected relativistic effects of nucleon motion in the 
deuteron are insignificant corrections (~ 0(1 — a) 2 )- Another potentially possible source 
of difference between nonrelativistic and light-cone descriptions is the nucleon spin rota- 
tion effect [j34f| . However it was demonstrated in Ref.p2| that in the case of perpendicular 
kinematics this difference is negligible once again. Note also, that relativistic effects are 
even less important for the FSI, since integrals in the FSI amplitude (see e.g. eq.flTB|)) 
are sensitive to rather smaller values of the deuteron internal momenta. 



5.2 Uncertainties in knowledge of the deuteron wave function 

Possible uncertainties in the calculation of the d(e, e'p)n processes with polarized and 
unpolarized targets are small and under control for the kinematics considered here be- 
cause nucleon momenta in the deuteron do not exceed 300 — 350MeV/ c and FSI terms 
are sensitive to smaller deuteron internal momenta. Another reason why our processes 
are less sensitive to the uncertainty of the deuteron wave function, is that we consider 
observables which are the ratios of experimental quantities in (nearly) similar kinematical 
conditions (see eqs. (|5~TD and (|54])). 
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5.3 Off shell effects 

The target nucleons are bound in the deuteron; the square of their four momentum is not 
the square of their mass. The nucleons are off-shell. One needs to estimate the influence 
of this effect to calculate the cross section of eq. ([□[). We estimated the uncertainties due 
to these effects by considering the Born amplitudes of models which account differently 
for off-shell effects in a ep |35], |22 . The differences are very small, less than one percent 



or so. 

One of the uncertainties originates from the deformation of the bound nucleon wave 



function. This deformation has been calculated in Ref.p6| for processes dominated by 
PLC. A similar effect arises within Skyrmion models of the two nucleon interaction [^7]]. 
The major effect is the suppression of the probability for PLC in bound nucleon due to 



the color screening phenomenon^]. The influence of this effect can be estimated by 



rescaling the deuteron wave function for a nucleon with momentum k by the factor: 

*(*)= fl + W 2 -^)-(l-f)"^^) , (56) 

where is the deuteron binding energy and AE(ph 0.6 GeV) is the parameter which 
characterize bound nucleon excitations in the deuteron. The Q 2 dependence accounts 
for the presumed dominance of PLC in bound nucleons for sufficiently large momentum 
transfer Q 2 > Q 2 ^ 2(GeV/c) 2 [^§. 

In Fig. 10 we present calculations of the ratio defined in eq. (|5l"l) (Fig. 10a) and the 
asymmetry defined by eq.(Q) (Fig. 10b) which include the influence of the PLC suppres- 
sion effect of Eq. (|5*B]) . Figures demonstrate that considered uncertainties are on the level 



of 5% and 10% in the case of unpolarized and polarized measurements respectively. 
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5.4 Semiclassical approximation for FSI in the momentum 
space - Feynman diagram approach 

When calculating FSI for the d(e,e'p)n reaction we substituted 0(— z) contribution in 
eq.(||) by \ (cf. Eq. ([])). This replacement allows to represent the FSI amplitude (Fig. lb) 
as the convolution integral of the deuteron wave function and the on-energy shell iViV 
scattering amplitude in the momentum space. The aim of this subsection is to substan- 
tiate our approximations. We start by analyzing the Feynman diagram, Fig. lb, which 
describes the FSI. Our interest is in the kinematics where the momentum of the nucleon- 

2 

spectator - p n is small: ^ -C 1. In this case it is reasonable to consider the nonrelativistic 
motion of the nucleon in the deuteron: % < 1, where p' n is the momentum of spectator 
in the intermediate state (p' n = p n — k, k is transferred momentum in the amplitude of 
FSI - see Figl.b). So it is legitimate to evaluate the loop integral by taking a residue over 
the energy of the spectator nucleon in the intermediate state - p' n0 - (For nonrelativistic 
motion of the nucleon this is the only pole in the lower part of the complex plane in the 
variable p' n0 .) Neglecting systematically all the terms ~ ^\ as compared to 1 we obtain: 



M'« _ (apt f^gP^M^ (57) 



<pVrlP> 2 1 A,-A, + k (27t 

Here the momentum q is chosen to be in the z axis direction, s = (p^ + q^) 2 and: 

Pnz = (x- l)m p=H2= = p nz - A = p nz - (E n - m) Md + q ° = (58) 

where p nz is the z-component of measured spectator momentum. The last term of eq. (|58|) 
accounts for real kinematics of the considered process. The amplitude f in the eq. ([57]) 
is normalized according to Imf NN (s,t = 0) = o to t and the deuteron wave function 
normalized as / ip 2 (p n )d 3 p n = 1. 

Eq.(|58|) reflects an important property of two-body high energy processes in which 
the variable = ko — k z (fig. lb) but not k z is small. This dynamics is easily accounted 



for in light-cone mechanics of deuteron [21 
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It is easy to demonstrate that eq.(|57|) is equivalent to the usual Glauber approximation 
if off-energy shell effects are neglected in f NN . Really, if we use the Fourier transform of 
the deuteron wave function as: 

1 



V»(z4) 



(27T 



(59) 



the integral in eq . (|57|) over p' nz differs from for z < only. This is the factor Q(—z) in 
the coordinate space representation. 

We shall neglect in the following analysis the dependence of / on p' nz ,p n z, since 
they are small for two body processes at high energies and will ignore off shell effects for 
small momenta of spectators. Then the integral over p' nz can be evaluated by deforming 
the contour integral over p' nz into the lower part of the complex plane. To calculate this 
integral we need to know analytic properties of the wave function of deuteron. For this we 
use the conventional parametrization of deuteron wave function calculated with Paris p0| 
and Bonn|TT| NN potentials: 



^(p) = E 



P 



(60) 



where J2j Cj 



3 r ' ' 3 

and s and d-waves differ by coefficients Cj. Substituting ip(p' n ) in 



eq. (|57D from eq.(|60D we obtain: 



M 



FSI 



(2tt 



< p\J e mJp > 



E 



d 2 p' n ± fd 2 P 



nz 



(2tt) 2 J (2tt) 



X 



/ 



NN 



(P'nz + ^PnX + m2 j)(P'nz ~ ^PnX + m ]) V ™ ~ P ™ + 16 
d'p'n. 



E 



(27T) 



NN 



C: 



(61) 



where p n = (p nz ,p' n ±)- Separating the real and imaginary parts inside of [...] and using 
eq.floTf) one obtains a factor | which we used in eq.(|) and additional term neglected in 
the above calculations: 



M 



FSI 



< p\ J e m^\p > 



fpL 

(2k) 



-f 



NN 



iPnz 



E 



C; 



2 2 

j 



Pn±+ m2 (pl + m2 ) 



(62) 
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where 



Pnz I 



± fNN 



(63) 



NN d2 P' n ± 
(2tt)2 



The size of (3 is a measure of the accuracy of replacing the Q(—z) function by the factor 
|. Eqs. (p^) , ( p5D and ( |5"3"D show that within the Glauber approximation the contribution 
from the factor f3 decreases with spectator longitudinal momentum p nz and is practically 
negligible for perpendicular kinematics. However in the case of the modified Glauber 
approximation for the three resonance model (sec. 3. 2. 3) the longitudinal momentum of 
intermediate state does not coincide with the longitudinal momentum of external de- 
tected particle. Due to this difference, the value of f3 in the three resonance model is 
determined at p' nz = jP nz — A m p , where A m p is defined according to eq.([45|). Therefore 
the contribution of (3 will be larger than within the conventional Glauber approximation. 
Note that contribution of A mj) in soft FSI are effectively included into the sum rule 



It follows from eq.(|62"D that the contribution of /3 factor to the interference of Born 
and FSI amplitude is further suppressed because the real part of f NN amplitude is small 
(~ 0.2 of the imaginary part) and because the Born term corresponding figl.a is real. 
The contribution of (3 factor to the double scattering term is ~ [3 2 (cf. eq. fl62|) ). 

In Fig. 11 we demonstrate the Q 2 dependence of the ratio defined by eq.(^l|) (Fig.lla) 
and the asymmetry defined by eq.(|54|) (Fig. lib) within the conventional Glauber ap- 
proximation and the modified Glauber approximation for the three resonance model. It 
follows from this calculation that including (3 within the three resonance model of modi- 
fied GA brings the predictions of this model and QDM approximation (where we included 
the 0(— z) factor explicitly in coordinate space) closer together. 
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5.5 Meson exchange currents (MEC) and A-isobar contribu- 
tions 

Estimates of contributions from meson exchange currents corresponding to diagrams 
similar to Fig.lc are rather controversial at large Q 2 since these terms are very sensitive 
to the assumed t-dependence of the meson-nucleon vertex form factors. These form 
factors are not obtained from theory; instead they are used as fitting parameters ||38||. On 



the other hand the restriction to x = ~ 1 and Q 2 > 1 (GeV/c) 2 strongly suppresses 
the sea quark content in nucleons. Including mesonic components is one way to treat the 
anti-quark content, so that a suppression of the sea can be considered as a hint of the 
suppression of MEC at large Q 2 and x > 1. This argument relies on the Bloom-Gilman 
duality [^] between the structure function of a nucleon at x — > 1 and contribution of 



resonances and on the fact that this hypothesis describes reasonably near threshold data 
for deep inelastic eN scattering. Besides, the diagram corresponding to Fig.lc contains 
an extra ^ factor as compared to the Born (Fig. la) and the FSI (Fig. lb) diagrams. 

The A-isobar contribution is expected to be small also in the kinematical range suit- 
able for studying CT effects, since the 7*iV — > A transition form factor decreases more 



rapidly with Q 2 than the N* transition form factors[40|. Another reason for the suppres- 
sion of the contribution of A's is that the AiV — > NN amplitude is predominantly real 
and decreases rapidly with energy (since it is dominated by pion exchange) whereas the 
FSI effects we study are determined by imaginary part of the soft rescattering amplitude. 

In fig. 12 we compare predictions of the Glauber approximations of this paper with 
results of calculations of the model of Ref. ||10||. [] Calculations has been performed in 
the kinematics Q 2 = 1 (GeV/c) 2 , qo ~ 400 — 500 MeV/c and spectator momentum 



p n = 400 MeV/c. Within model of Ref. Jl0| the contribution of meson currents and 



isobars in the kinematics where NN rescattering dominates is small (~ 6% for MEC 
and ~ 4% for isobar contribution). Comparision of our approach and that of Ref. |10 1 



4 We are thankful to W. Leidemann for making available his calculations within the approach of 
Ref. flOf] for the kinematics discussed in this paper. 
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hints that Glauber approximation is applicable for Q 2 as small as 1 (GeV/c) 2 . Note that 
calculation according to Ref . [|K| can be considered as the upper limit for the contribution 



of meson currents. Theoretical analyses of current data (cf. Ref. |[41||) and calculations of 



Ref. indicate that the contributions of meson exchange currents have been strongly 



overestimated previously. 

5.6 Charge exchange contribution and spin dependence of the 
elementary amplitude 

In this analysis we have neglected contributions of terms in which the large longitudinal 
momentum is first transfered to the neutron which then converts to a proton via a 
charge exchange process. Since the charge exchange reaction np — > pn is dominated by 
the pion exchange its amplitude is predominantly real. Therefore it mostly contributes 
to the double interaction term, not to the screening term. We can estimate its relative 
contribution as 

^chargeexchange ^ &en,elastic d(T V P j ffy (64) 
^double @ep,elastic 



Using the data from Ref. [4~3|, we estimate that this correction is of the order of 9% 
for Q 2 ~ 2GeV 2 and it decreases rapidly with Q 2 , approximately as Q~ 4 . For the case 
of leading neutron production, the relative contribution of the charge exchange is larger 
by a factor o 2 ja 2 en < 5. Note that we expect transparency effects to cut down the 
charge exchange just as they reduce elastic scattering. This is the chiral transparency 



effect discussed in Ref. f24jl ; so this charge exchange effect does not affect our conclusions 
concerning the sensitivity of the d(e, e'p)n reaction to color coherent effects. 

Note also that we neglected another effect of dependence of the elastic pn amplitude 
on the parallel or anti-parallel nature of the helicities of the colliding nucleons. If we 
consider, for certainty, magnetic transitions, the rescattered nucleons will predominantly 
be in the antiparallel (parallel) state for = 1(0). Using the current data on the 
o^ijm) difference |[44|| , we estimate that this effect leads to a correction to the rescattering 
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amplitude which is < 3% for Q 2 ~ 2(GeV/c) 2 and < 1% for Q 2 ~ A(GeV/c) 2 . This effect 
is somewhat more important for Ad calculations. We will consider the effects discussed 
here in more detail elsewhere. 

5.7 Factorization approximation 

In the DWIA form of cross section (eq.(0)) we assume that the electron- nucleon cross 
section is not changed by the distortion of kinematics due to final state interaction. This 
commonly used assumption for the distorted wave impulse approximation allows us to 
factorize the cross section in the form of eq.(|l2|). However computing the unfactorized 
form of cross section for the deuteron is quite straightforward and the authors will discuss 
it elsewhere. For small nucleon momenta within deuteron this effect is small anyway. 

6 Summary 

We have demonstrated that selecting perpendicular kinematics for the spectator: p n ~ 
200 — 400 MeV/ c, x ~ 1 significantly increases the effect of struck nucleon rescattering 
off the spectator nucleon. Thus varying p n allows the selection of smaller than aver- 
age internucleon distances in the FSI. If p n is small the final state interactions cause 
a screening effect, and the effect of including the color transparency is to increase the 
computed cross section over the value obtained using the Glauber approximation. But 
for larger values of p n , the conventional final state interactions enhance the cross sec- 
tion. In this case, the influence of color transparency is to suppress the computed cross 
sections. Based on these observations, we suggest measuring the ratio of d(e, e'p)n cross 
sections, where in one case FSI is dominated by screening effect and in the second case 
by the double scattering. This ratio is more sensitive to the effects of color transparency. 
Our calculations of this ratio, using several CT models, predict 20 — 40% CT effects for 
Q 2 ~ 4 - 10 (GeV/c) 2 . 

High sensitivity of the d(e, e'p)n cross section to the final state interaction of suffi- 
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ciently energetic struck nucleoli is one the of preconditions for looking for color trans- 
parency phenomena. This condition is valid in the case of high energy electrodisinte- 
gration of a polarized deuteron with production of a spectator under the perpendicular 
kinematics discussed above. We suggest measuring the tensor asymmetry in d(e, e'p)n 
reaction, where conventional Glauber approximation predicts a huge influence of FSI. 
The dominance of deuteron d-state wave in the polarized density matrix restricts the 
essential internucleon distances to ~ 1.5 fm, and suppresses the PLC evolution to a 
normal size nucleon. The CT models we consider give predictions of 50 — 100% effects 
at Q 2 ~ 4 - 10 (GeV/c) 2 . 

To estimate the possibility for the unambiguous investigation of CT effects with 
deuteron target we considered a number of theoretical issues which might affect the 
reliability of the theoretical calculation of the d(e, e'p)n reaction. The analysis shows 
that under our perpendicular kinematic conditions one can confine the theoretical un- 
certainties to the level of ~ 5% and ~ 10% for unpolarized and polarized measurements 
respectively. 

Comparison with models explicitly accounting for meson currents and isobars shows 
that in the perpendicular kinematics preferable for searching for CT effects the deuteron 
can be considered as a clean detector with a low rate of noise. 
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Figure Captions 



Figure 1. Graphs for d(e,e'p)n scattering, (a) - Born (PWIA) approximation; (b) - final 
state interaction contribution, where the broken line accounts for all NN multiple 
scatterings; (c) - meson exchange current contribution. 

Figure 2. The dependence of transparency - T (eq.(|49|)) on the spectator neutron angles, for 
different values of neutron momenta: (a) p n = 100 MeV/c, (b) p n = 200 MeV/c, 
(c) p n — 300 MeV/c and (d) p n = 400 MeV/c. The solid-straight line is the Born 
approximation, dashed line - contribution of interference between Born and FSI 
amplitude, dash-dotted - effects of screening - is the sum of Born and interference 
terms, dotted line - contribution of the double scattering term, solid-line - is the 
sum of all terms in the cross section. Q 2 = 6 (GeV/c) 2 . 



Figure 3. Q 2 dependence of transparency T defined according to eq. (|49|) at x = 1. Solid line 

- GA prediction, dashed line - CT prediction within QDM approximation, with 
AM 2 = 0.7 GeV 2 and dotted line PWIA prediction. Curves labelled "I" are the 
calculations at p n = 0, "II" - at p n = 100 MeV/c, "III" - at p n = 200 MeV/c "IV" 

- at p n = 300 MeV/c, and "V" - at p n = 400 MeV/c. 

Figure 4. Q 2 dependence of transparency T defined according to eq.(^) normalized on GA 
calculations: dash-dotted line GA approximation, dotted line - QDM prediction at 
p n = 200 MeV/ c, solid line - QDM prediction at p n = 300 MeV j c and dashed line 

- QDM prediction at p n = 400 MeV jc. The curves labelled by boxes correspond to 
Q 2 dependence of QDM predictions for R(Q 2 ,p n i,Pn2) defined according to eq. (|5T|) , 
normalized on GA calculations. Solid line with "boxes" - p n \ = 300 MeV/c, p n 2 = 
200 MeV/c and dashed line with "boxes" - p nl = 400 MeV/c, p n2 = 200 MeV/c. 
QDM predictions calculated with AM 2 = 0.7 GeV 2 . x = 1. 

Figure 5. Q 2 dependence of transparency T at p n = 200 MeV jc - (a) and p n = 300 MeV jc 

- (b), for x = 1. The solid-line - prediction within the GA, dotted line - QDM 



35 



approximation with AM 2 = 0.7 GeV 2 , dashed line three resonance model within 
modified GA, with parameters M N * = 1.4 GeV, M N „ = 1.8 GeV, e = 0.17, 
J N ' N = 1.0 and ,f iV *' iV = 3.1. Solid and dashed lines labelled by boxes are 
calculations within the GA and three resonance model of CT within Green function 
method and zero-range NN scattering approximation. The parameters for three 
resonance model are same as above. 

Figure 6. Transverse momentum dependence of unpolarized density matrices - Pd{k, k) = 
u{k) 2 + w{k) 2 (eq. flTBl) ) - dashed line and polarized density matrices - P2o(k,k) 
(eq.(|53"D) - solid line, at p z n = 0. 

Figure 7. The (9 n ,p n ) dependence od asymmetry defined by Eq.(^). (a) - Born (PWIA) 
approximation, (b) - Glauber approximation. Q 2 = 6 (GeV/c) 2 . 



Figure 8. Q 2 dependence of asymmetry defined according to eq.(j54|), at x = 1. Solid line 



Glauber approximation, dashed line - QDM of color transparency, dash-dotted line 
- the three resonance model of color transparency within the modified GA. Straight- 
dotted line - Born approximation, (a) - p n = 250 MeV/c, (b) - p n = 300 MeV/c, 
(c) - p n = 350 MeV/c and (d) - p n = 400 MeV/c. QDM and three resonance model 
parameters are same as in Fig.5. 

Figure 9. Q 2 dependence the ratio R(Q 2 ,p n \-,Vn2) defined according to eq. (fTC]) , at p n \ = 
300 MeV/c and p n 2 = 200 MeV/c - (a) and asymmetry defined according to 
eq. (|54l) at p n = 300 MeV/c - (b). x = 1. Solid line - Glauber approximation, 
dashed line - QDM approximation with AM 2 = 0.7 GeV 2 , dotted line - AM 2 = 
0.9 GeV 2 and dash-dotted line - AM 2 = 1.1 GeV 2 

Figure 10. The theoretical uncertainties in determining of Q 2 dependence of transparency 
R(Q 2 ,Pni = 300, p n2 = 200) - (a) and asymmetry A d (p n = 300), at z = 1. 
Solid-line - glauber approximations, dashed - lines QDM approximations. Upper, 
down solid, dashed lines at Q 2 = 1 (GeV/c) 2 - are the GA and QDM predictions 
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with Paris and Bonn wave functions respectively. The curves inside corresponds 
to predictions with Paris light cone wave function, color screening effects (which 
according to eq . (|56|) appears at Q 2 > 2 (GeV/c) 2 ) and predictions with different 
approximation for a en (eq. (|T2D . QDM calculated with AM 2 = 0.7 GeV 2 . 



Figure 11. Q 2 dependence the ratio R(Q 2 ,p n i,Pn2) defined according to eq.([5ll), at p n \ = 
300 MeV/c and p n2 = 200 MeV/c - (a) and asymmetry A d defined according to 
eq.(|54|) at p n = 300 MeV/c -(b). x = 1. Solid line - Glauber approximation, 
dashed line - three resonance model of CT, within modified Glauber approxima- 
tion, doted line - QDM approximation of CT, solid line with "boxes" - Glauber 
approximation with (3 factor considered in eq.(|62|) and dashed line with "boxes" 
- three resonance model of CT, within modified Glauber approximation and with 
(3 factor corresponding to higher mass excitation in the intermediate state. QDM 
and three resonance model parameters are same as in Fig.5. 

Figure 12. The 9 n dependence of ratio R. Curve marked by "GA" shows the ratio of d(e, e'p)n 
cross section calculated within Glauber approximation, without factor of f3 (eq. (|63|) 
to the PWIA cross section, " GA"' - same as above, with Glauber cross section cal- 
culated with factor (3, "NM" - ratio of cross section with multiple scattering and 
Siegert operator to the Born cross section, which also accounts for the direct pro- 
ductions of slow neutrons, " NM+MEC" - " NM" calculation with meson exchange 
currents, and "NM+MEC+IC" -"NM" calculation with meson exchange currents 



and isobar contribution. The last three curves are those of the calculation of Ref . fllP 1 
for p n = 400 MeV/c and Q 2 = 1 (GeV/c) 2 . 
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